Mon. Not. R. Astron. Soc. 000, [THT5] (2012) Printed 9 March 2012 



(MN WFeX. style file v2.2) 



The UV peak in Active Galactic Nuclei : a false continuum 
from blurred reflection ? 



O 
(N 



A. Lawrence^ 

Revised version : March 2nd 2012 

^ Institute for Astronomy , SUPA (Scottish Universities Physics Alliance), University of Edinburgh, 
Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK 



Accepted 2012 March 8th. Received 2012 March 2nd ; in original form 2011 July 18 



OO 

o 



(N 
OO 

o 



ABSTRACT 

I summarise and analyse key problems with observations of the UV bump in Active 
Galactic Nuclei (AGN), and especially the accretion disc interpretation - the temper- 
ature problem, the ionisation problem, the timescale problem, and the co-ordination 
problem - and suggest that all these problems can be solved if, in addition to the ac- 
cretion disc, there is a population of cold, thick clouds at approximately 30i?s which 
reprocess the intrinsic continuum. Exploring cloud parameter space, I find that clouds 
with density n ~ lO^^cm"^ and column Nh > 4 x lO^^cm"^ reflect most of the in- 
trinsic continuum, but convert a substantial fraction of the EUV luminosity into lines, 
dominated by Ly/3 and Hell Lya. When velocity-blurred, this makes a false continuum 
peak at '^llOOA which fits the observed SED well, but turns back up in the FUV to 
make a hard EUV SED, as required by ionisation models. I argue that the observed UV 
variability is dominated by this component of fixed shape, possibly due to changes of 
covering factor. The amount of mass required is small, so it is not necessary to disrupt 
the disc, but only to make an unstable and inhomogeneous atmosphere. The proposed 
clouds may be related to those suggested by several X-ray phenomena (X-ray reflection 
components, high-velocity outflows, Compton thick partial covering) but are not the 
same, leading to a picture with a wide range of inhomogeneous structures at different 
radii. 
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1 INTRODUCTION 

The most prominent feature in the Spectral Energy Distri- 
bution (SED) of Active Galactic Nuclei (AGN) is the so- 
called "Big Blue Bump" (BBB ) . This broad feat ure peaks 
in the UV I Sanders et al. I ll989l : lElvis et al.|[l99i ). roughly 
consistent with the expectation of thermal emission from gas 
accreting onto a black hole mass with a mass of 1O^~®M0. 
The BBB is however broader than a single black body, as 
one would expect if accretion takes place i n a disc with 
a range of temperature s at different radii ('Shicldi 1 19781 : 
iMalkan fc Sargent! [iQsl) . In more detail, modellers have 
found it hard to reproduce the shape of the SED, especially 
near the UV peak, and through the sof t X-ray region (e.g. 
iBlaes et al.ll200ll : iHaroCorzo et al.ll2007^ . 

Some observations provide generic support for the idea 
of a disc-like structure in the centre of AGN. (i) The strong 
Fe line and hard X-ray flattening seen in the X-ray spectra 
of AGN are often interpreted as being due to reflection by 
"cold" material th at covers a solid ang le of ~ 2n as seen by 
the X-ray source IPounds et al.lll990l ). (ii) In a handful of 



quasars, optical spectra in polarised light, removing contam- 
inating emission, show a broad dip possibly corresponding 
to the Balmer ed ge absorption expected from an accr etion 
disc atmosphere I Kishimoto. Antonucci fc Blaesll2003l ). Ex- 
tending the polarised continuum into the near-IR reveals the 
classic long wavelengt h u^^^ spectrum expect ed from simple 
accretion disc models l|Kishimoto et al1l2008l ) . (iii) For radio 
loud objects, and especially superluminal objects, where we 
can have some indication of the orientation of the sources, 
AGN show velocity width a nd and surface brightness corre- 
lations with viewing an g le llBrothertonlll996l : iRokaki et al.l 
l2003l : IJarvis fc McLurd ^0061 ) . (iv) At optical wavelengths, 
over a very broad range of luminosities AGN show colour- 
luminosity trends consistent with th e generic expec tation 
from multi-temperature black bodies (|Lawrencell2005h . 

It is notable that all the above successes concern op- 
tical or near-IR wavelength emission. In the UV however 
a number of observations disagree badly with at least the 
simplest accretion disc models. These issues can be charac- 
terised as the temperature problem, the ionisation problem, 
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the timescale problem, the co-ordination problem, and the 
size problem. Some of these problems, along with other key 
iss ues such as polar isation predictions, have been discussed 
bv lAntonuccil (|2002l ) and references therein. In the next sec- 
tion, I look at the current state of each of these key problems 
in turn. In the following section, I consider an idea that has 
often been seen as a rival to accretion disc models - emission 
from dense clouds. In fact, a combination of accretion disc 
and dense inner clouds may be exactly what is needed to 
solve the problems of accretion disc models, as large veloc- 
ity blurring can produce a localised false peak in the SED. 
Finally, I briefly discuss implications of this idea, and how 
it can be tested. 



2 PROBLEMS WITH THE BIG BLUE BUMP 

Here I summarise a number of problems with our current 
understanding of the Big Blue Bump. 



2.1 The temperature problem 

AGN seem to be cooler than they ought to be. The SEDs 
of AGN seem to show a universal near-UV shape, reach- 
ing a maximum in uSv at a wavelength of ar ound llOOA, 
more or regardless of luminosity or redshift |Zheng et al.l 
19971: iTelfer et al.ll2002l : IScott et al.lliooi : IShang et al.il2005l : 
Binette et al. 200^~Such a peak suggests a characteris- 
tic temperature of T~ 30,000K, and indeed this is similar 



to the inner temperatur e of accretion disc fits ever since 
iMalkan fc Sargent! l| 19821 ). However, for any thermal model, 
the characteristic temperature should be roughly 



Tch 



where Le ~ L/LEdd is the luminosity in units of the Ed- 
dington luminosity, Mg is the mass of the central object in 
units of l(f Mq, and 7?5 is radius in units of ^Rs, where Rs 
is the Schwarzschild radius. Here -Rs is taken to be the ra- 
dius of an optically thick spherical surface from which black- 
body radiation at temperature T emerges. Other estimates, 
such as the radiation-density equivalent temperature, will 
be similar, and other possible characteristic temperatures, 
such as the Compton heating temperature, wi l l be h otter, 
and rather model dependent (|Ferland fc ReesI Il988l ). The 
maximum temperature of a simple accretion disc extend- 
ing to the last stable orbit is substantially hotter than the 
naive temperature above. Discs in which one imposes the 
zero-stress condition at the inner boundary can be cooler in 
the middle and hotter towards the outer radii, but it is far 
from clear whet her this applies t o black hole accretion discs 
(see for example iFrank. King fc R ainc 2002, p. 90). That de- 
tailed accretion disc models do not overcome this problem 
is clear fo r example f r om th e attempts to fit the SED of 
3C273 bv iBlaes et"al] Hool). Smaller black holes are hot- 
ter, but objects radiating at much less than the Eddington 
limit are cooler, so for local Seyferts, as opposed to luminous 
quasars, we might expect some diversity; but these also seem 
to show the universal knee. 

Fig. [1] shows the SED of the nearby quasar 3C 273, 
plotted in linear units, compared to a blackbody of tem- 
perature T~32,000K which is a g ood first approximatio n to 
the observed SED. According to IPeterson et al.l (|2004h the 



mass of the black h ole i n 3C 273 is 8.87x10^ Mq . (Note that 
iKaspi et al] (|2000l ) and lPahani fc Tiirlej l|2005f ) find values 
roughly a factor of two either side of this value.) Integrat- 
ing over the SED shown in Fig. [TJ and extrapolating the 
FUV with a = —1.8, the total luminosity is approximately 
Lbs ~4x10''^W, so that Lb ~0.356. If this luminosity were 
radiated from 5Rs, the expected black body temperature 
would be T~76,000K. This is clearly discrepant with the 
observed SED, as also shown in FiglU 

One could fix the temperature problem by making AGN 
radiate well below the Eddington limit, but only by mak- 
ing their black holes almost two orders of magnitude more 
massive. Alternatively one could require that the radiation 
emerges from a larger radius - of the order 50-Rs . Of course 
most of the energy generation has to take place interior to 
this radius, or once again an unreasonably large mass is re- 
quired. One is led therefore to the possibility that interior 
energy is absorbed and reprocessed, and that the peak we 
see corresponds to the reprocessed rather than the original 
emission. 

The FU V SED shortward o f the llOOA peak i s mor e 
controversial. IZheng et all ([l993) and iTelfer et al] l|2002h . 
by compiling a composite SED using HST data spanning 
a large range of redshifts, find that the FUV falls steeply 
all the way to 300A as oc v" with a ~ —1.8. Con- 
structing such a high-z composite requires careful correction 
for Lya for est absorption by the intervening IGM. On the 
other hand, IScott et al] (|2004h use FUSE data for relatively 
low redshift AGN to construct a composite reaching 650A 
and find a mean slope of a ~-0.6. They find considerable 
FUV diversity, and a tendency for high-lumi nosity AGN to 
be stee p, and low- luminos ity A GN to be flat. Binette et alj 
( 20051) . IShang et al.l (|2005h . andEaroCorio et al. I (|2007l ) ex- 
amine the FUV SEDs of individual AGN, up to z~l, where 
the IGM correction is still s mall, and confl r m the diversity of 
FUV shapes. In particular, iBinette et al.l (|2005l ) stress that 
while some quasars continue to decline in vL^ out to 600A 
others turn back up at short wavelengths. Fig. [T] shows ex- 
amples of these two cases taken from their paper. 

Another way of altering the expected emi ssion is by ab- 
sorpti on. Various possibilities are reviewed bv lBinette et alJ 
(|2008l ). The observed knee is close to the Lyman edge, but 
the shape is not like an edge, and it is definitely at some- 
what longer wavelength. It is also possible that dust red- 
dening has altered the shape - de-reddening the 3C273 SED 
by Av ~ 0.35 over and above the line of sight Galactic ex- 
tinction of Av ~ 0.1 ma kes the match with the expected 
temperature much better. iBinette et al.l (|2005l ) explore the 
possibility that crystalline dust can produce a broad absorp- 
tion trough matching that seen in PG1008-I-1319 and similar 
objects. Overall, absorption models remain plausible, but I 
do not consider them further in this paper. 



2.2 The ionisation problem 

The ionisation problem is a corollary of the temperature 
problem : an SED falling in vL^ terms does not seem capable 
of producing the broad emission lines observed in AGN. This 
is especially true i f the high-redshif t com posites derived by 
IZheng et al.l (ll997D and lTelfer et"all (I2002D are correct; these 
show a very soft FUV continuum, with a ~ —1.8, whereas 
the observed emission lines in the optical and near-UV seem 
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Figure 1. Comparisons between various observed SEDs and simple theoretical curves. Note that both axes are li near, which helps t o 
emphasise how strongly the SED is dominated by the NUV peak. The dotted lines show the SED of 3C 273 taken from lKriss et al\ \l99f\ ). 
kindly supplied by G.Kriss, using nea rly simultaneous data from various instruments. The datap oints have been de-re ddened assuming 
E(B-V)=0.032 (see \Kriss et al\\l99!\) . The circles show continuum from t wo quasars analysed bv lBinette et al\ t200A) . estimated from 
their figures. The dashed line labelled "T02" shows the composite SED of \Telfer et al\ \200^) , represented as a power law with index 
-1.8. The blackbody curve with T=32,000K is chosen as an illustration of how a temperature in this region dominates the SED. The 
blackbody curve wth T=76,000K represents the temperature that would naively apply g iven the observed lum inosity and mass of 3C 273 
(see text). The dotted line labelled "K97" represents the ionising continuum used bv lKorista et al\ tl99'i) . The normalisations of all 
components are arbitrary; they are scaled for comparison purposes. 



to indicate a much harder continuum, with a ~ —0.5. 
Fig. [1] contrasts the observed SED with the ionisin g contin- 
uum u sed in the comprehensive BLR study of Kori sta et al.l 
l|l997l ) - a power law with index a = —0.5 and an exponential 
cut off at A;r=44eV - which does a good job of reproducing 
the observed ratios of broad emission lines in A GN. The ioni- 
sation problem has been known for many years (Netzer 1985; 
CoU in-Souffrin 1986; Dumorit, C olhii- Souffrin & Nazarova 
199^1 : iKorista. Ferland fc Baldwinil 19971 '). Netzer's (1985) ar- 
gument concerned the ratio of the total BLR emission to 
Lya. Roughly speaking, if Lya tells us the continuum near 
the Lyman edge, then extrapolating with a steep slope does 
not provide enough energy to produce all the BLR emission. 
Netzer argued that the problem could be solved if the BLR 
(bu t not the conti n uum) is reddened hy Av ^ 0.25-0.5, 
but ICollin-SouffrinI | 198a ) shows that the observed contin- 
uum then_^jb_to_£roduce_thc reddening-corrected line ra- 
tios. iKorista. Ferland fc Bald win (199^) concentrate on the 
Hell lines, especially HeII1640, and sho w that a steeply ex - 
trapolated continuum in the manner of IZheng et al] (| 19971 ') 
fails to produce the observed Hell 1640 equivalent width by 
at least a factor of 5, and probably considerably more. Ko- 
rista et al consider several possible solutions. One is that 
the UV-EUV bump is double peaked. This is the possibility 
that I pursue in this paper. 

Fig. [T] illustrates the ionisat ion problem , contr asting 
the observed SE P of 3C 273, t he iTelfer et al.l (|2002l ') com- 
posite, and the IKorista et al.l (|l997 ) ionising continuum. 



Interestingly, although the observed SED seems inconsis- 
tent with both the required ionising continuum, and the 
naively expected temperature, the ionising continuum and 
the naively expected temperature are roughly consistent 
with each other. Furthermore the FUV upturn shown by 
PG1008-I-1319 is roughly consistent with the slope expected 
from the naively expected temperature, suggesting that this 
very hot emission is present, but somehow masked in the 
NUV region. 



2.3 The timescale problem 

Relatively few AGN have well sampled UV light curves, but 
those that do show fairly well defined ups and down s. The 
situation is summarised bv lCollier fc PetersonI (|200ll ). who 
derive UV structure functions for four AGN (NGC 3783, 
NGC 7469, NGC 5548, and F9) and optical structure func- 
tions for six other AGN. They show characteristic timescales 
(defined by the fiattening of the structure function) of be- 
tween 5 and 94 days. Collier and Peterson show evidence 
that this characteristic timescale is proportional to black 
hole mass, which for these UV-studied local AGN is of the 
order lO'^ - IO^Mq. This is roughly confirmed by the fact 
that much more lumi nous quasars seem to show timescales 
of the order of vears (iGiyeon et ai1ll999l: iTurler et al.|[l999 



Trevese. Kron fc Bunond l200ll : iTrevese fc Vagnettil bood : 
Hawkins, 20071 

A timescale of the order of tens of days is in- 
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consistent with simple quasi-steady accretion disc mod- 
els, as the viscous (radial drift) times c ale is far longer, 
as was first recognis e d by lAlloin et al.l l|l985h . Following 
iFrank. King fc Raind (|2002l ), Ch 8, and scaling to the time 
expected for a radius appropriate to a standard thin disc ra- 
diating at a characteristic NUV temperature of T=30,000K, 
the viscous timescale is 



ivis 



12.6 yrs x L] 



0.1 Mo.i 



where a and /i are the usual viscosity and efficiency param- 
eters, scaled to their expected values. Acoustic modes also 
do not work. For T=30,000K, assuming the gas is largely 
ionised, the sound speed is ~ 22 km s~^, so that the sound 
crossing timescale is t sound = 12.6 yrs xiiaoMg. On the other 
hand, the light crossing timescale, tught ~ 8.3 hrs x iZaoAfg 
is much too short, although it may be relevant to the co- 
ordination of variability at different wavelengths, as dis- 
cussed below. 

There are two timescales that may work. The first is 
some kind of dynamical timescale. The free fall timescale 

3/2 

from 30Rs is tff ~ 5.9 days x R^^ Mg. This is slightly on 
the fast side, but very much in the right ball park, sug- 
gesting that accretion flow may be dynamically unstable. 
The second possibility is some kind of cooling timescale, 
which will for example determine the time in which inter- 
nally stored energy could be released as radiation in a flare. 
This is much more model dependent, depending on the total 
mass involve d, its temperature, and w hether it is optically 
thick or thin. lCoUier fc PetersonI (J200l) state that for a stan- 
dard accretion disc model the thermal or cooling timescale is 
tthermai ~ 27.5 days X R^Q^ Msag\ where a is the standard 
disc viscosity parameter. 

If either a dynamical timescale or a thermal timescale 
is the correct answer, the variability cannot be modelled by 
a quasi-steady accretion disc - the physical structure of the 
disc is changing faster than the gravitational energy is being 
dissipated. However this problem could be removed if what 
we see is reprocessed emission - only the reprocessor has to 
change fast, not the disc itself. In section 14.31 we will look 
at whether dense reprocessing clouds can provide the right 
timescale. 



2.4 The co-ordination problem 

The next serious problem for accretion disc models is 
that variations at different UV-optical wavelengths occur 
in phase - i.e. the peaks and troughs line up, so that 
any delay or smearing is small compared to the variabil- 
ity timescale itself, which is very hard to achieve if the ra- 
diation at different wavelengths come from different parts 
of the disc. The best evidence comes from a handful of 
local Seyfert galaxies with well sampled UV light curves 
fNG C 5548 : iciavel et al.lll99ll: iKorista et al.|[l995l : NGC 



4151 :lEdelson et al.lll996l : ICrenshaw et al.ll996l: N GC 7469 : 
iKriss et al.ll2000l : NGC 3516 i lEdelson et al.ll200 g'). For more 
luminous quasars, there are few UV datasets, but fairly well 
sampled long term optical light curves presented by Giveon 
et al. (1999) and Hawkins (2003) show changes at differ- 
ent wavelengths that are also aligned to much less than 
the typical variability timescale. The lack of wavelength- 
dependent delays led to the suggestion that the inner ac- 
cretion disc is heated at least in part by the central X-ray 



source (rather than locally by viscous torques), so that the 
UV variations are actually a repr ocessed version of varia- 
tions in the central X-ray source l|Krolik et al.lll99ll ). The 
expected delay between different wavelengths is then the 
light travel time between different regions, which is of the 
order hours. Comparing UV wavelengths, several authors 
have c laimed a measurement of such a short de l ay for NGC 
7469 jWanders et al.l Il997l : ICollier et al.l 1 19981 : iKriss et all 
|2000| ). but for other AGN (NGC 5548, 3516, 3783, 4151, 
and Fairall 9), attemp ts have only resulted i n limits of the 
order 0.1 to 0.3 days jC rcnshaw ct al. ■1996"; E delson et al.l 
1 19961 . I2OOOI : iPeterson eTal. T998 ). Korista fc GoadI (l200ll ) 
have claimed that the NGC 7469 effect is an artefact 
due to contamination of the continuum by diffuse emis- 
sion from the BLR. Comparing different optical wave- 
lengths, a number of authors have shown marginally signif- 
icant delays of the order of days in several different AGN 
llCollier et all I2001I: lOknvanskii et~al1 l2003l: fSergeev et all 



I2OO5I : iDoroshenko et al.ll2005l : iBreedt et al.ll2009h . 

A prediction of the "disc heated by X-ray source" pic- 
ture is that the UV and optical emission should lag the X- 
ray emission. However, simultaneous monitoring campaigns 
show that any such delay is less than 0.15 day s in both NGC 
4151 (|Edelson et al.l 1 19961 ) and NGC 3516 l|Edelson et al.l 
i2000l ). and in the case of NGC 7469, if anything t he UV 
leads the X-ray emission by ~ 4 days. [Sreedt et al] l|2009h . 
from 5 years of monitoring MKN 79, find that overall. X-ray 
and optical emission are well correlated, with an optical vs 
X-ray lag of less than ~2 days, but also note differences. The 
X-ray variability has more short timescale power, as might 
be expected, but when smeared with the expected disc trans- 
fer function it fails to reproduce the observed optical light 
curve (see their Fig. 7) - the predicted light curve matches 
some but not all of the observed flaring on timescales of 
tens of days, and does not produce the long term trends on 
ti mescales of hundreds of days. S imilar problems were noted 
bv lKazanas fc NavakshinI (|200ll ). 

A further problem is that the X-ray luminosity in AGN 
is much smaller than the the UV luminosity. For relatively 
low luminosity AGN, the total X-ray luminosity is sensi- 
tive to how the observed 2-10 keV spectrum is extrapolated, 
so that X-ray lu minosity may in fact be comparable to the 
UV luminosity (|Chiang fc Blaej 12003 ) . However, this can- 
not be a universal explanation, as the X-ray /UV ratio is an 
order of magnitude s maller for the most luminous quasars 
(Stra teva et al.l [20051 ). Overall there is clearly a close con- 
nection between UV-optical emission and X-ray emission, 
but the relation is a complicated one. One possibility is 
that, rather than the X-ray emission being primary and the 
UV a reprocessed version, the real primary emission is the 
FUV peak at ~ 300A which is predicted by simple accre- 
tion disc theory and photo-ionisation calculation. We do not 
have FUV light curves, but one AGN, NGC 5548, has been 
monito red the other side of th e possible peak, at ~ 80 A with 
EUVE. [Marshall et all ([l993) showed that over ten days the 
variations at 76A and 1350A are well correlated, with any 
lag less than 0.25 days, but that the amplitude of varia- 
tions was several times larger in the EUV. In a later EUVE 
observation made s imultaneously with RXTE and ASCA, 
I Chiang et al. I l|2000l ) showed that the 76A emission leads the 
X-rays by 0.4 days. There is a case to be made then that both 
the UV and X-ray emission are substantially reprocessed. 
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2.5 The size problem 

Sizes deduced from the fluctuations of gravitationaUy lensed 
quasars seem to be discrepant with expectation. The light 
from each lensed component can suffer additional micro- 
lensing by stars, or possibly by dark matter substructure. 
The expected optical depth is close to unity, so that each 
component has ongoing irregular fluctuations that are un- 
correlated with each other, and can therefore be distin- 
guished from intrinsic variability co mmon to all compo - 
nents. (See for example the review of IWambsgansj \200d ) 
The magnification depends on source size, as the expected 
size, while small compared to the Einstein radius of the 
lenses, is not point like - i.e. a larger source has more 
washed out fluctuations. A clear qualitative result is that 
component-to-component anomalies are significantly bigger 
in X-rays than in the optical-UV, showing that the X-ray 
source is more com pact, by a factor of several (|Dai et al.l 
I2OI0I : IChartas et al,,2010. ) and possibly t hat the hard X-ra y 
size is smaller than the soft X-ray size (|Chen et al.ll2012l ). 
Likewise, broad emission lines show little or no fluctua- 
tion, and therefore the BLR size is substantiall y bigger than 
the continuum source size (|Lewis et al.lll99"^ ). Within the 
optical-UV r ange, fluctuations are l arger at shorter wave - 
lengths, (e.g. lEigenbrod et al.|[20oi : iMosquera et al.ll201ll ). 
strongly supporting the general idea of a distributed source 
with a colour gradient. The slope is plausibly consistent with 
the radial temperature gradient expected from an accretion 
disk, but with a rather large error. This is an area that will 
hopefully improve strongly in coming years. 

As well as these striking relative results, lensed quasar 
monitoring can give us information on the absolute source 
size. This is somewhat model dependent as it requires un- 
derstanding the macroscopic lensing, the mass function of 
the microlenses, and the relative transverse motions. Early 
result s claimed either consistency with expecte d source 
sizes l|Wambsganss. Schneider fc Paczvnskil Il990l ) or that 
obser ved sources sizes were sma ller than expected accretion 
discs ijRauch fc Blandford|[l99ll ). More recent and substan- 
tial studies have consistently concluded that observed sizes 
are a factor of several larger than expected from simple 
thin a ccretion disk m o dels o f the same objects (IPoolev et al.l 
I2OO7I : [Morgan et all I2OI0I : I Jimenez- Vicente et all 120121 ) . 
while correlating clearl y with the blac k hole mass estimated 
from line widths (Mor gan et al.lfioiol ). Suggested solutions 
to this problem could include much lower accretion effi - 
ciency, or a flatter temperature profil e |Morgan et al.ll2010l ) 
or an inhomogeneous accretion disc (iDexter fc Agoll|201lf ). 
The central reprocessing region discussed later in this paper 
is another possible solution to this problem. 



2.6 Amplitude effect and mixing model 

Related to the co-ordination problem is the fact that the 
amplitude of variability for both local A GN and luminous 
quasars is strongly wavelength depend ent (|Cutri et al.ll 1985 1: 



MJ: 

IClavel et aLlll99ll : IWilhite et aj1l2005l'). This i s mos t clearly 
demonstrated in Fig. 13 of IWilhite et all |2005| ). which 



shows the ratio of the mean difference spectrum to the mean 
spectrum for 2500 SDSS quasars. This variability index in- 
creases steeply at wavelengths shorter than 25QQA. A generic 
way to explain the amplitude and co-ordination effects at 



the same time is by mixing a constant red component with 
a variable blue component (|Lawrence 2005 '). However, only 
some such models work. Fig. |5] shows two deliberately sim- 
ple toy mod e ls app lied to the NGC 5548 campaign data of 
IClavel et al.l l|l99ll ). The constant component is modelled as 
a power law with slope a and normalisation NpL, and the 
variable component as a blackbody at a single temperature 
T with normalisation Nbb- 

In the first model, the temperature T of the black 
body was fixed, and flux-ratio values were calculated for 
a sequence of values of the black-body normalisation Nbb, 
to produce predicted tracks in two diagrams - the ratios 
F1337/F1825 and F1337/F267Q vs the flux F1337, where 
F1337 refers to the flux at 1337A . The other parameters 
{NpL, OL, and T) were adjusted to produce a good simultane- 
ous (eyeball) fit between these tracks and the data. The best 
fit temperature for the variable black body was T=27,500K, 
with the power law of the fixed component having slope a 
= -2.0. 

In the second model, the emitting area of the black body 
was held fixed, and the temperature T varied to produce 
the predicted tracks. This of course produces variation in 
both colour and normalisation. As with the first model, the 
free parameters {NpL, a, and black body emitting area) 
were adjusted to try to achieve a simultaneous fit to the 
two flux ratios, but no flt could be found. As an illustration. 
Fig. [2] shows parameters adjusted to flt the run of the ratio 
F1337/F1825, which results in a prediction for F1337/F2670 
which is a long way off. 

Of course both models are physically unrealistic - the 
variable component is unlikely to be a single blackbody, the 
red component may not be well represented by a power law, 
and the variable component may vary in both area and tem- 
perature. Nonetheless, the difference between the two simple 
extremes is striking. It strongly suggests that the way to ex- 
plain the variability of NGC 5548 is by the variability of a 
component with fixed shape. 



3 REPROCESSING BY DENSE INNER 
CLOUDS 

3.1 Introduction 

The temperature problem (i.e. the premature turnover of the 
SED) and the variability behaviour of AGN suggest that 
the observed UV continuum represents re-processed emis- 
sion. The earlier idea that the UV represents reprocessed 
X-rays seems untenable because there is not enough X-ray 
luminosity, and because in at least some cases the X-ray 
emission lags the UV emission. The primary radiation could 
instead be the unseen peak EUV emission - the inner disc 
heating the outer disc. Against this idea is the observation 
of steeply declining FUV continuum. In support of the idea 
is the fact that emission lines suggest that there is indeed a 
hard EUV continuum, and the fact that at least some high 
redshift quasars show a turn up shortwar d of 800 A. Per- 
haps, as suggested by iKorista et al.l |l993), the continuum 
is double peaked. 

However, there remains the puzzle that the wave- 
length of the vLi, peak (~ llOOA) seems to be universal 
(|Shang et al.ll2005l ). even if the continuum a few hundred 
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Figure 2. Mixing models for the colour- a mplitude tren d in th e 
variability of NGC 5548. The data are from lClavel et al\ ^199i) . 
The open circles show the flux ratio F1337/F2670, and the filled 
circles show F1337/F1825, where F1337 refers to Fx at 1337 A. 
The model tracks show the effect of a fixed power law and a vari- 
able black body. In the left hand pane, the temperature of the black 
body is fixed at T=27,500K, with only its normalisation varying. 
In the right hand pane, the emitting area of the black body is fixed 
and its temperature varied. 



A shortward of this is not (Scott et al.l [ 20041 : iBinette et alJ 
I2OO5I : iHaroCorzo et al ] [2003)7 Furthermore, within individ- 
ual objects, it seems that multi-wavelength changes may be 
explained by the variability of a component of fixed shape. 
Where a spectral shape is relatively constant across a large 
range of conditions, this is suggestive of a solution depen- 
dent on atomic physics rather than macroscopic physics. 
Perhaps what we see as the UV Bump is actually made 
of emission lines. The features that we recognise as the 
"broad emission lines" come characteristically from a radius 
of ~1000Rs. Clouds at say ~30Rs would produce emission 
lines 6 times broader, which would be seen as a false con- 
tinuum. A similar logic has led some previous workers to 
consider optically thin (primarily free-free) emission from 
hot (10^~'')K gas as an alternative to the accre tion disc 
ijFerland. Korista fc Petersonlll990l : lBarvainislll993t ). but my 
aim here is to discuss the role of optically thick clouds in 
conjunction with the accretion disc. In the next subsection, 
I discuss some general issues concerning dense inner clouds. 
Following this, I show model results for emission from such 
clouds, how they could produce a false local continuum peak 
at 1100 A , and discuss the prospect for more detailed mod- 
els. In section |4] I discuss whether we expect such clouds to 
exist, and how this relates to other issues such as the X-ray 
spectra of AGN. 



was either non-thermal (with power law slope -1) or a very 
hot X-ray plasma, with the issue being whether "cold" ther- 
mal material could survive within the non-thermal plasma, 
and whether the reprocessed thermal emission could be a 
viable alternative to an accretion disc for explaining the 
optical-UV emission in AGN. These ideas seem to have 
gone out of fashion because it has become widely ac- 
cepted that an accretion disc is the primary mode of en- 
ergy generation. An overlapping series of papers has con- 
sidered the effect of such cold clouds on the X-ray spec- 
trum, especially the Fe line and the Compton reflection 
hump llLightman fc Whit3 Il988l: ISiyron fc Tsurutal Il993l: 



Karas et al.l I2OO0I : iMalzac fc Celottil I2OO2I : [Merloni et all 
20061 ) ■ and also the possibility that variab le partial cover - 
ing explains X-ray variability ( Abrassart fc Czernvl l2000l l . 
However this has been a minority industry, as most workers 
have followed [Pounds et all l[l990l ) in assuming that X-ray 
reflection is from the surface of the accretion disc. 

My approach in this paper differs from the above pa- 
pers in two respects. Firstly, rather than trying to produce 
an alternative to the accretion disc, I will examine how cool 
dense clouds can modify the accretion disc SED. Secondly, 
I will concentrate on a somewhat different region of param- 
eter space. Most of the above papers consider very dense 
clouds in the very central regions (3Rs<R<10Rs) where 
most of the energy is generated, arguing that they can be 
confined by the strong magnetic fields that may be present. 
[Kuncic. Blackman fc Reeg|[ 19961 ) argued that such clouds 
must be very small and dense (hydrogen density n^lO'^* 
cm~'^Rj^M^^) if they are to cool in less than the free-fall 
time. The radiation from such clouds will be close to black 
body at around lO^K, as shown by detailed calculations in 
nearly all the above papers, and as such would be hard to 
distinguish from other forms of cold material, such as an 
accretion disc, and indeed this emission would suffer from 
many of the same problems previously discussed in section[51 
However, at somewhat larger radial distances, say 30Rs, the 
density required to cool in a reasonable time is much less, 
n^lO'^^. Such clouds will radiate mostly in atomic transi- 
tions, but they will still be moving fast enough that the 
Doppler blurred emission will produce a false continuum. 
If the accretion disc itself is the source of such material, 
then there is no confinement issue, and the starting tem- 
perature of the clouds (^lO^K) will be similar to the ex- 
pected photo-ionisation equilibrium temperature when ex- 
posed to the ~10^K radiation from the inner region. I will 
consider the origin of the clouds a little more closely in sec- 
tion |4] For now, I will assume that moderately dense clouds 
exist at R^SORg, and consid er the consequences. I note 
that ICzernv fc Dumont[ ([l998l ) have also examined n-^lO^'^ 
clouds, but they consider clouds at rather higher ionisation 
parameter, and assume a non-thermal primary continuum. 



3.2 Cold dense clouds 



3.3 Emission from cold dense clouds 



Cold dense clouds in AGN were a minor fashion for 
just ov er a decade (iGuilbert fc Ree3 [1988: Fe rland fc Rees 
19881 : [Lightman fc White! [ 19881: [Celotti. Fabian fc Rees 



I992I: [Sivron fc Tsuruta[ [ 19931: [Colhn-Souffrin et al.[ [ 19961: 



Kuncic. Blac kman fc R ees 'l99 6|; Kunci c. Celotti fc ReesI 
1997. : .Czcrnv fc Dumont.,1998, : iCelotti fc Rees 199^ ). In aU 
these papers, it was assumed that the primary energy source 



I examined the properties of clouds ionised by a hot central 
source at a 3C 273-like luminosity, using ver sion C08.00 of 
the w ell known "Cloudy" software package ([Ferland et al.l 
[l998h . Except where noted below the ionising source was 
a blackbody with temperature r=100,000K and (total) lu- 
minosity L=4xl0^^erg s~^, and the distance of the clouds 
from the source was J?=10^®cm, corresponding to 30Rs for 
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M_rf = lO^M0. The calculations were performed assuming an 
open geometry. A variety of densities and column densities 
were examined, as described in the subsections below. The 
temperature of 100,000K could be seen as an approximation 
to the disc emission interior to the clouds, or if perhaps the 
disc is disrupted, could be thermal emission from the very 
dense interior clouds discussed in section 13.21 In many of 
the cases described below, I tried varying the temperature 
over the range 75,000 to 200,000, at fixed total luminos- 
ity, and found no qualitative difference. In a few cases, I 
also tried a continuum which added a two component X- 
ray continuum to the 100,000K black body, based on the 
data in Kriss et al (1999). I used a soft X-ray power law 
with slope a=-1.5, a low energy break at 7.5 Rydberg, and 
luminosity L=5xl0'*^erg s~^; and a hard X-ray power law 
with Qf=-0.6 and luminosity L=1.25x 10*^erg s~^. The main 
effect of such an X-ray continuum was to turn a back-end 
neutral zone into an extended partially ionised zone. The 
effect on the reflected spectrum described below was signifi- 
cant but not large - around 20%. Finally, the choice of 30Rg 
is guided by the logic described above in section [321 but I 
also considered some clouds at other distances, as described 
below. 

In the subsections below I refer to logn and log(A''H) 
where the hydrogen number density n is in units of cm^"^ 
and the hydrogen column density Nh is in units of cm^^. 

3. 3. 1 BLR-like clouds, log n ~ 10 

Clouds with density logn=10, such as are often assumed 
in BLR models, can be considered for column densities 
log Nh <24.5. Above this Nh value, the cloud size is larger 
than the whole region. Throughout this column density 
range, the clouds are fully ionised in H and He. At low 
columns, log Nh < 23, the clouds are essentially transpar- 
ent. Above this range, the electron scattering depth becomes 
significant, and an increasingly large fraction of the contin- 
uum is reflected back. The total emission, which is what 
one would see if equal numbers of front faces and rear faces 
are visible, is not very different from the input continuum. 
There are significant emission lines, but not strong enough 
to produce the effect we are looking for here. Clouds at this 
density may of course exist in the inner region, and pro- 
duce observable absorption and emission lines, but will not 
produce a gross modification of the SED. 

3.3.2 Dense clouds, logn ~ 12 

Clouds with density logn=12 are much more promising. At 
low column densities, they are ionised, with the results being 
fairly similar to lower density clouds and with little net effect 
on the SED. However at log A^j^ =24.55 and above the rear 
part of the cloud becomes neutral. (With X-ray continuum 
added, this rear zone is partially ionised, but the transition is 
in the same place.) As the total column increases above this 
value, the ionised column stays the same and the neutral 
column simply increases. The ionised column corresponds 
to the point where the electron scattering optical depth is 
2.3. Varying the density from logn=ll to 13 changes the 
transition optical depth a little, from 4.0 to 0.9. 

The emission from the cloud is more or less the same 



for all columns log Nh >24.6, and is very asymmetric. Seen 
from the far side, the cloud is essentially opaque at all wave- 
lengths. However, the cloud is a strong reflection source in 
both continuum and lines. The resulting spectrum is shown 
in Fig. |3l The electron scattering depth of the ionised zone 
is of order unity, so that a large fraction of the continuum 
is reflected back. However around 8% of the incident EUV 
continuum is reprocessed into emission lines, the strongest 
lines being He II Lya at 304A and HI Ly^ at 1025A. The 
summed luminosity of the lines in the lOOOA region is ~30% 
of the incident luminosity at <lRyd. This, together with 
the marked drop shortward of the Lyman edge, is sufficient 
when velocity-blurred to make a substantial but localised 
false continuum, as discussed in section [3741 

Over the range 22< log A'^^f <24.5 the luminosity in Ly/3 
increases by more than two orders of magnitude, as shown 
in FigU while the ratio Ly/3/LyQ increases from ~ 0.1 to 
~ 8, and the ratio Ly/3/He II Lya increases from ~ 0.1 to ~ 
1.4. From log A''£f =24.5 to the maximum plausible column of 
log Nh='27, the luminosity in Ly/3 is constant, and the ratios 
Ly/3/Lya and Ly/3/He II Lya are constant at 8.51 and 1.38 
respectively. If we fix the column density at log Nh='25 and 
vary the temperature of the ionising source from T=75,000 
to T=200,000, the ratio Ly/3/LyQ changes little (from 8.71 
to 8.32), but the ratio Ly/3/He II Lya changes from 4.27 to 
0.45. This last sensitivity may give a way to test models of 
this kind. 

Clouds with logn=ll give a similar result. Fig.[5]shows 
how the luminosity in various lines changes as the density 
is varied, for logAH=25 and T=100,00K. The Ly/? effect 
is significant over a range of density, but strongly peaked 
around logn ~12. 

Similar clouds much closer in give a somewhat different 
result. At 5Rs, a cloud with logn=12 and log Ah=25 still 
has a neutral back end, and so is opaque when seen from the 
rear, but the larger ionised column results in larger electron 
scattering depth, t^s = 5.3, so that the reflected continuum 
is little different from the incident continuum, and the re- 
flected lines are a order of magnitude smaller than the 30Rs 
case. Further out, similar clouds at 100-Rs have res=0.7 and 
give a result similar to the 30-Rs case, with Ly/3 in fact sev- 
eral times stronger, although of course these lines will be 
distincly less broad than those at 30-Rs. By 300-Rs, line re- 
flection is still strong, but now dominated by Lya rather 
than Ly/3. 

3.3.3 Very dense clouds, logn ~ 14 

What about clouds with even higher density, say log n > 14? 
This strays into the region which Cloudy was not designed 
for, and so issues warnings, but the rough results are clear. 
At low columns, such clouds are ionised, but at logA''i^ >23 
they are mostly neutral with a thin ionised skin. There is 
little continuum reflection; the reflection spectrum is domi- 
nated by HI Lyman continuum emission, with He I Lyman 
continuum and He II Lya also strong, and significant HI 
Balmer continuum. The relative strength of these compo- 
nents seems to be roughly constant at log Nh >24. I do 
not explore these densities in more detail, because the re- 
sults from Cloudy may be unreliable and also because they 
do not seem to give the 1100 A peaking that observations 
require. 
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Figure 4. Dependence of reflected line luminosity on col- 
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are calculated on the basis of a J^-k covering factor. 
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Figure 5. Dependence of reflected line luminosity on density. 
In all cases the other parameters are \ogNH=25, T=100,000K, 
log Ltatal=46. 6 andlog R=16. The luminosities are calculated on 
the basis of a J^n covering factor. 



At even higher densities, clouds will of course start to 
look like black bodies, and emit at a temperatures in the 
range 10*~^K, depending on radial distance, as discussed 
by several previous authors (see discussion in section 

3.4 Prospect of dense cloud reflection models 

Adding velocity blurring to the reflected spectrum of a sin- 
gle cloud can produce a false continuum of roughly the right 
kind. What about an ensemble of clouds surrounding the ac- 
cretion disc - will the emission be dominated by the intrinsic 
continuum or the reflected spectrum ? This will depend on 



both the covering factor of the clouds, and the amount of 
blurring. 1 blurred the spectrum shown in Fig.[3]with a Gaus- 
sian of FWHM 75,000 km 'sT^ , roughly what one might ex- 
pect for clouds at 30-Rs, depending on the kinematic model. 
The reflected quasi-continuum vSv at 1100 A is then roughly 
twice the incident continuum at that wavelength. If the cov- 
ering factor is C then the amount of leaked-through intrin- 
sic continuum is proportional to 1 — C and the fraction of 
the total due to reflection, ignoring self-occultation effects, 
is then approximately 2C/(1 -|- C). For the reflection spec- 
trum to dominate therefore requires C > 1/3. The fraction 
changes quite quickly around this value, with C=0.1 giving 
18%, and C=0.5 giving 67%. Larger covering factors require 
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self-occultation to be taken into account. There is there- 
fore only a small range of covering factors which produces a 
strong reflection spectrum without self-occultation. 

However self-occultation, with large covering factor, 
may be exactly what we need. Note that the reflected spec- 
trum from a single cloud, like the incident continuum, is 
dominated by emission at several Rydbergs. Clouds which 
see only other clouds, rather than the disc surface, will also 
produce a spectrum like that of Fig. [S] The cloud spectrum 
bounces from one cloud to another and emerges at the edge 
of the system of clouds at a kind of quasi-photosphere. The 
central disc will be completely occulted, so the SED that 
we see is composed of the pure cloud spectrum plus the un- 
occulted outer part of the disc. Fig. [6] shows a toy model of 
this kind compared to the combined data for 3C 273 and 
PG 10084-1319. The model shows a blurred version of the 
Fig. [3] reflection spectrum, added to a u^^^ power law with 
an exponential cut-off at 0.4 Rydbergs, meant as an approx- 
imation to an accretion disc truncated at 30Rs- The rela- 
tive strengths of the two components are roughly correct, 
but note that this is not a full or self-consistent model. It 
should be seen as a plausibility demonstration rather than a 
real model. Significant extra work, and decisions on model 
assumptions, is needed before it makes sense to attempt a 
proper fit. 

(i) There will very likely be a wide range of densities and 
column densities at different radii. To some extent, the SED 
modification may naturally be dominated by clouds in a par- 
ticular range - by comparison with these just-right clouds, 
clouds of lower density or column may be transparent, and 
clouds that are closer in may primarily reflect an unaltered 
continum. However, a prediction of the actual range of cloud 
parameters could make a big difference. 

(ii) The shape of the SED modification will be sensitive 
to the radial distribution of velocity and cloud density. In 
this paper, I make basic arguments about the likely density 
and velocity regime concerned. Going beyond this needs a 
physically motivated model which makes a prediction of the 
velocity and density distributions. 

(iii) A realistic model for the intrinsic continuum is also 
needed, both for the interior region, which provides the con- 
tinuum which ionises the clouds, and the exterior region, 
which produces the optical SED. The behaviour of the in- 
ner accretion disc is clearly crucial. 

(iv) Geometrical considerations could also make a big 
difference. For example, it may be that clouds lifted out of 
the disc in a narrow range of radii will form a kind of funnel, 
somewhat akin to the model of Elvis (2000) but much closer 
in. In this case, the continuum that we see will be inclination 
dependent. 

(v) To fit the whole optical-UV-FUV SED, we may need 
to take into account a more continous system of clouds, 
between my suggested ~30R5 region and the traditional 
~1000Rs BLR region. 



clouds, the origin of UV variability, and the implications for 
AGN spectra. 

4.1 Size of clouds and overall cloud system 

Clouds with density n^f^lO^^ cm~^ and column N_ff=10^^ 
cm""' will have thickness D~10^"^cm. If spherical, they will 
have individual masses of the order 1O~''M0, but of course it 
is quite possible that the material is the form of much more 
extended sheets or filaments. The total mass of the system 
of clouds can be estimated by considering a spherical sheet 
at radial distance R with covering fraction C, which gives 

TUs = 8.3M0 X CMgRioN25 

where as usual Mg is black hole mass in units of 1O^M0, 
R30 is the radial distance in units of 30Rs, and N25 is in 
units of 10^^ atoms cm~^. How does this compare with the 
mass likely to be involved in the accretion flow onto the 
black hole ? Let us first consider simple uniform spherical 
accretion. Assuming steady flow in free-fall the density at R 
is 

n = 2.2 X 10** cm"^ x fio\Mg^LERio 

where Le = L/LEdd, and /j, is the usual radiative efficiency 
scaled to an expected value of 0.1. The accreting mass inte- 
rior to R is 

rua = I.IMq x fio lMlLERlij^ 

The mass available is therefore an order of magnitude less 
than the required mass of clouds, and a very large degree of 
clumping would be needed to make clouds of the necessary 
density, so it seems that a simple spherical accretion model 
does not work. 

There are two ways to increase the mean density and 
mass available. The first is for accretion to be inefficient 
{^J. <^ 1), as in ADAF models, where most of the energy gen- 
erated by infall is advected into the black ho le. Potentially a 
clum py ADAF model could be constructed (|Celotti fc ReesI 
[l99i), but I do not consider this further here as it seems 
unlikely to be relevant to luminous quasars. 

The second way to increase the mass available is to 
slow down the infall and so increase density by a fac- 
tor K = Vff/vinf. In standard accretion disc models 
K ^ {R/H)^ where H is the characteristic scale 

height of the disc and a is the usual viscosity parameter 
(|Frank. King fc Raineir2002l . Chapters 5 and 8). The density 
is further increased by the vertical compression factor R/H. 
For an AGN disc, assuming a ~ 0.1 and R/H ^ 10, then the 
mass available within R is at least rria ^ 10'^ MqM^^ LeR%^ 
and the mean density in the disc at 7? is at least n ~ 
10 Mg Le cm~ i?3Q , with each of these quantities likely 
to be orders of magnitude larger. In conclusion there is a 
plentiful supply of material in the disc. We look at disc mod- 
els a little more closely in the next section. 



4 DISCUSSION 

Do we expect that such dense thick clouds will exist in the 
inner regions of quasars ? Here I briefiy discuss the mass 
of the possible cloud system, its origin, the lifetime of the 



4.2 Possible disc origin of clouds 

Detailed disc models (eg iHubenv et al.l I2OOC ) confirm the 
above rough calculation. From the figures in Hubenv et al.l 
1 estimate that for a disc radiating at 15% of Ed- 
dinston around a black hole of mass IO^Mq, the mass in- 
terior to 30Rs is 7x10'^Mq, so only a very small fraction 
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of the disc mass is needed. The density varies greatly, horn 
~ lO'^^ cm""^ in the midplane to ^ 10^" cm"'' at the disc 
surface. Likewise the total column density is of the order 
10^* cm~'^ or more, depending on radius, and thinning out 
towards the surface. Production of clouds then does not re- 
quire complete disruption of the disc, but just a clumping 
of the upper parts. It has been argued that in the inner re- 
gions discs become unstable to clumping (Krolik 1998) but 
in fact only a modest degree of clumping would be required 
to produce the necessary clouds. 

To provide the reflection medium we are considering, 
the clouds also need to be lifted out of the disc, and to 
be cool enough to be not already completely ionised, i.e. 
a few times lO^K. Discs are expected to have a signifl- 
cant vertical temperature gradient, so again this means that 
the upper portions would be the source of cloud material, 
rather than the disruption of the whole disc. In disc wind 
models (see eg iProgal |2007| and references therein) , uplift 
and gas temperature are closely linked. Radiation pressure 
comes from absorption of resonance lines as well as elec- 
tron scattering, so that even well below the nominal Ed- 
dington limit gas can be lifted above the disc surface, where 
it is then hit by the central UV source and driven out. In 
the very inner disc, it is generally argued that the mate- 
rial is too highly ionised for this effect to work, so that 
winds are launched at ^ 100_ Rg, and only for high mass 
and h i gh luminosity systems (eg Shlosman. Vitello &i Shavivl 
1 19851 : iMurrav et al] Il995l : iProga. Stone fc KallmanI l200d ). 
However. iProgal (|2005*) shows that in the inner regions a 
"failed wind" or "puffed up disc" can be produced, which 
may be exactly the kind of structure we need. Note that 



local radiation force is very sensitive to clumping. If clumps 
can form within the upper disc with the kind of parame- 
ters we have considered here, log n=12 and log Nn—'ib, they 
will be optically thick to the Lyman continuum, so that the 
"force multiplier" will be extremely large. Such clumps are 
then very likely to be ejected from the disc when formed. 
A reasonable guess is that at large radii (> 100-Rs) gas can 
be uplifted without clumping, and so form a smooth wind, 
whereas at very small radii, (< 10-Rs) the material is too 
dense and hot to either feel uplift forces from heavy element 
absorption lines, or to make the cool thick clouds needed to 
get Lyman uplift. The result may be that dense clouds are 
formed and uplifted at a characteristic radius around SORg. 
More detailed calculations are needed to see if this is indeed 
the case. 

Alternatively, clouds maybe uplifte d by magneto- 
centri fu gal effects, as in the models of iKonigl fc Kartjg 
(1994). IdeKool fc BegelmanI (|l995l ). and lE verettl l|2005f ). 
This idea is usually considered at much larger radii, but 
may also be relevant in the inner regions. 

4.3 Cloud timescales and variability 

Do we expect such clouds to survive ? Unless confined (eg 
by an external medium or magnetic field) they would dis- 
rupt on the sound crossing timescale. Assuming that the gas 
is mostly ionised and at a temperature T30 = T/SOOOO this 
gives tdiarupt = 52.0 daysx N25nl2T^f^^'^ , which is very sim- 
ilar to the dynamical timescale at 30-Rs. Because there is a 
plentiful supply of material, we do not require the clouds to 
be long lived. The cloud sound crossing timescale will also 
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be (very roughly) the timescale on which clumps could grow 
inside the disc, so we may expect that clouds are repeatedly 
forming and being destroyed. 

How quickly would such clouds be expelled ? If the 
clouds originate in the disc, they will initially be moving 
at orbital velocity, but because they are optically thick 
to Lyman continuum and electron scattering, they will be 
driven radially outwards with a greatly super-Eddington 
force. Given that they have mass per unit area Nnmn, 
the time to reach escape velocity will be Tesc = 124 days x 

3/2 1 

seems therefore that clouds will be de- 
stroyed faster than they are expelled, and so we do not nec- 
essarily expect a wind to form - more of a broiling cloudy 
atmosphere. 

The covering factor of clouds may well vary significantly 
on the cloud formation/destruction timescale, leading to a 
natural explanation of UV variability. The intrinsic UV con- 
tinuum, and the optical continuum, could vary quite slowly; 
it is just the reflected continuum that varies rapidly. The 
timescale is about right for the observed UV variability, and 
the fixed shape of the reflection spectrum is consistent with 
the strongly wavelength dependent variability, and the lack 
of delays between wavelengths. There could still be hours- 
timescale light travel time delays, just as in the disc re- 
processing model. A possible objection to the idea that the 
intrinsic disc emission varies slowly is that the BLR emis- 
sion lines, which track the observed UV continuum with a 
time lag, must surely be responding to the EUV continuum, 
which then must also be varying rapidly. However the dense 
inner clouds also reflect EUV continuum; the BLR clouds 
will see a mixture of slowly varying intrinsic EUV continuum 
and rapidly varying reflected EUV continuum. Whether this 
works quantitatively should be explored. 

4.4 Effect on X-ray spectrum and variability 

Much has been written about absorption and reflection ef- 
fects in X-ray spectra (see the review by [Turner fc Millei] 
l2009l 'l. Ultimately a self consistent model for both UV and 
X-ray emission is needed, which is beyond the scope of this 
paper. I restrict myself to a few general points. 

(i) In the picture emerging in this paper, the inner ac- 
cretion disc is assumed to be present, so blurred ionised disc 
reflection will certainly occur. However, the clouds we have 
proposed will also produce blurred X-ray reflection effects, 
albeit at somewhat lower velocity blurring than the inner 
disc effects. If the inner disc is occulted, the cloud reflection 
effects may dominate. The clouds will also produce opaque 
partial covering, which can e asily lead to an appa rent reflec- 
tion efficiency greater than l. lMerloni et all (|200d 'l have con- 
sidered inhomogeneous X-ray reflection models, using dense 
clouds with no central accretion disc, but embedded in a very 
hot plasma. Their clouds are of similar column to those dis- 
cussed here, but at higher density, and closer in. They find 
that soft excess and hard X-ray humps naturally go together. 

(ii) Our clouds are close r and thicker than those de- 
duced bv lRisaliti eFail l|2007D and related papers to explain 
rapid absorption changes in many AGN. They may be more 
closely related to to the material producing Compton-thick 
partial covering phenomena revealed by BeppoS AX and 
Suzaku at E >10keV (Dadina fc Cappi, .2004: Reeves et al.l 
I2OO9I : [Turner et al.ll2009i r although it has been claimed that 



such objects may also be explained with disc reflection 
iFabian et al.l [2005'). From a stud y of 165 AGN with IN- 
TEGRAL data. iRicci et al.l l|201ll ') have proposed that the 
large reflection component seen in apparently mildly ab- 
sorbed Type 2 AGN may be explained by Compton-thick 
clumps partially covering the source. At middling cover- 
ing factors such opaque partial covering produces a striking 
spectral signature. At large, but still not complete, covering 
factors, the leaked-through spectrum may look quite nor- 
mal, and if there is additional thinner absorbing material at 
larger radii, the spectrum may easily be mistaken for be- 
ing Compton-thin when in fact it is largely Compton-thick. 
A consequence is that the observed X-ray luminosity may 
be systematically wrong. The effect on the luminosity func- 
tion for obscured and unobscured AGN could produce an 
art iflcial dependence of obs curation on luminosity, as noted 
bv lLawrence &l Elvid l|2010l ) and explored in more detail by 
Mayo and Lawrence (in preparation). 

(iii) The clouds we have discussed may also be re- 
late d to the high velocity absorbers seen in highly ionised 
Fe JPounds et al.ll2003l : IChartas. Brandt fc Gallagheill2003l : 



I Reeves et al.l 



■ .2009!^ to be very common 

i Tombesi et al.ll2010l ). and have typical velocities of 0.1-0. 2c. 
It is generally assumed that they represent outflows, but 
there have also been claims of inflows a t similar velocity 
(iDadina et al. I I2OO5I : iTombesi et ail 120101 1. If material is in 
rotational or turbulent motion at such velocities, this is quite 
consistent with the kind of bl urring we have ass umed in this 
paper. In at least one case (Risali ti et"ai1l2005l 'l the outflow 
velocity was seen to change substantially from one obser- 
vation to another, suggestive of the kind of instability we 
have suggested for our UV-reflecting clouds. The columns 
deduced for the high velocity X-ray absorbers are large but 
somewhat less than we have been considering here, log Nh = 
23 — 24 as opposed to log Nh > 25, but we note that partial 
covering will dilute the strength of the absorption lines. Our 
UV-reflecting clouds may also be related to the fast moving 



the ubiqitous soft excesses seen in AGN ( 


Gierlinski & Done! 


I2OO4I: iMiddleton. Done & Gierlinskill2007 


) 



(iv) In section l4^ I suggested that variations in covering 
factor could explain the UV variability. Likewise, covering 
factor changes co uld explain the obse r ved X -ray variabil- 
ity, as argued by lAbrassart fc Czernvl l|200Cf ) . Several pa- 
pers have argued that absorption changes can explain the 
com plex changes seen in the spect ra of individual objects 
(e.g. iMiller. Turner fc Reevea 20091 and refe rences therein; 
see also the review bv lTurner fc Milleij|2009l '). However, UV 
and X-ray variability cannot be explained by exactly the 
same reprocessing clouds. As discussed in section 12.41 X- 
rays and UV track each other well over medium-timescales, 
but X-rays also show extra short timescale variability. This 
is naturally explained if UV reprocessing is dominated by 
thick clouds at ~30Rs, whereas X-rays are affected both by 
these clouds and by closer-in clouds as well. However, it is 
not obvious why X-ray variations sometimes seem to lag the 
UV variations. 

Overall, the clouds required to explain Nh changes, X- 
ray variability and spectral changes, thick partial covering, 
soft excesses, and high velocity outflows are all likely to be 
somewhat related to our proposed UV-reflection clouds, but 
are not the same. A distribution of cloud properties is re- 
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quired, and it is not yet clear if a single model can reproduce 
everything consistently. 

4.5 Apparent UV source size 

Following equation 2 of lMorgan et alj l|201(]| ). the character- 
istic UV size of a thin accretion disc around a black hole 
with M = W^Mq accreting at the Eddington limit with ef- 
ficiency 77 = 0.1 should be ~ 10-Rs. Ifowever we have found 
that reprocessing clouds could be producing a false UV con- 
tinuum over a considerable radial scale, 10-100 Rs - energy 
produced in the central region emerges on a larger scale. 
Potentially therefore reprocessing clouds can solve the size 
problem implied by gravitational microlensing studies (see 
section [2.5^ . Whether this works quantitatively is not yet 
clear - the effective size will depend on the density and cov- 
ering factor of the clouds with radius, which obviously re- 
quires a physical model. Such a model of the physical distri- 
bution of the clouds is also needed to predict the chromatic 
behaviour of microlensing fluctuations. 

4.6 Next steps 

In this paper I have concentrated on the properties of clouds 
that could plausibly explain the worrying observed features 
of the UV bump and its variability, and demonstrating the 
plausibility of such clouds as an explanation. What is needed 
is a more detailed specific model. This is rather a challenge; 
to go beyond the plausibility demonstrated here one would 
need to specify the distribution of clouds in density, col- 
umn, and radius, and to specify their velocity field. This is 
too large a volume of parameter space to blindly explore as 
a grid. It really needs a prediction from a physical model. 
However, there are some general features that could be ex- 
plored first. 

(i) There should be some testable dependence on black 
hole mass and accretion rate, as these properties together 
predict the hardness of the intrinsic continuum. As noted in 
section 13.3.21 the relative strengths of Ly/3 and Hell Lya, 
and hence the 1100 A bump and the EUV peak, will depend 
on the temperature of the intrinsic continuum. This may 
produce measurable effects on the UV-FUV SED, and on 
the observed BLR line ratios. 

(ii) Another possible mass dependence comes from 
cloud size. The expected cloud size is ~ lO^^cm, but with 
a likely order of magnitude range. For quasar-size black 
holes, M ~ W^Mq, this is a small fraction of the source 
size, so that any one cloud has a small covering factor, and 
any changes are global ones pertaining to the ensemble of 
clouds. For lower-luminosity Seyferts, with Al ~ lO^'^Af©, 
the movement and physical changes of individual clouds may 
be significant, producing qualitatively different changes. 

(iii) Although I have argued that physically motivated 
cloud models are needed next, an obvious step is models with 
a small number of cloud zones, to see whether UV-reflection 
effects, X-ray reflection effects. X-ray absorption effects, and 
microlensing chromatic effects, can all be explained together 
without the various cloud populations destructively interfer- 
ing as it were. This could also lead to insights which would 
help the design of simultaneous UV / X-ray monitoring cam- 
paigns. 



5 CONCLUSIONS 

The problems with the accretion disc interpretation of the 
Big Blue Bump - the universal knee and FUV falloff, the 
ionisation paradox, the timescale issue, the co-ordination 
and amplitude-colour problems, and the microlensing size 
discrepancy - can all potentially be explained if EUV emis- 
sion from the inner disc is reprocessed by dense clouds at 
a few tens times Rs- A variety of clouds may be present, 
but the ones that produce the desired effect have densities 
in the range 10" < n < 10^3 cm , and column densities 
Nh > ■ix lO^"*. Such clouds are opaque and reflect most of 
the incident continuum, but recycle a substantial fraction of 
the EUV energy as emission lines, dominated by Hell Lya 
and HI Ly/3. Blurring this emission produces a localised false 
peak at llOOA, followed by a turn back upwards in the FUV. 
Producing such clouds requires only a small fraction of the 
mass available in the accretion disc, and does not require the 
accretion disc to be completely disrupted, and may be most 
likely to happen at a few tens of Rs . Such clouds will feel a 
very strong outward radiation force, but they may be formed 
and destroyed on a faster timescale than they are ejected, 
so that there will not necessarily be a continous wind, but 
rather an unstably variable cloudy atmosphere. Variations 
in covering factor may explain the observed UV variability. 

Such clouds may be related to those deduced from vari- 
ous X-ray observations, but are not the same. There is likely 
to be a large range of cloud properties at various different 
radii in an AGN, with those in different parts of parame- 
ter space producing different effects. It remains to be seen 
whether a physically realistic model can explain all the ob- 
servations simultaneously and consistently. 

In the 1980s, there was confused argument over whether 
the "blue b ump" was really due to an accretion disc as 
proposed by iMalkan fc Sargent (Il982l ) - some researchers 
claimed that the turn-up seen in optical spectra was 
actually due to Balmer continuum and Fell emission 
lines on top of an underl ying power-law (|Grandil Il982l : 
I Wills. Netzer fc Willsl 1 19851 ). The solution was to take a 
wider (IR-optical-UV- X- ray) perspective and see that both 
phenomena were present. Hi VIS fc Lawrence! (Il985l ) first used 
the term "little blue bump" to distinguish the atomic fea- 
ture from the "big blue bump" . If it is correct that the 1 lOOA 
peak is due to blurred atomic features, we may have to re- 
name this feature the "medium blue bump", with the true 
"big blue bump" residing at 300A. 
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